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Abstract — In this paper, we have a presented the cubic spline 
interpolation method for zooming in an image. The cubic spline 
algorithm is programmed in MATLAB as m-file. Firstly, 
according to the zooming intensity, sufficient zeroes are inserted 
between each pixel and then those zeroes are replaced by the 
respective values generated by the m-file. Finally a zoomed in 
image in constructed. 


f 

Si(x) if xi<x<x2 
Ete(x) if X2<X<X3 


Si(x)= < 


Index Terms — Digital image processing, zoom in, cubic spline 
interpolation method, MATLAB. 


Sn-l(x) if Xn-l<X<Xn 

V 


I. Introduction 

The field of digital image processing refers to processing 
digital images using digital computers. A digital image is 
comprised of finite number of values each of which has a 
location and a value. Generally speaking these values are 
referred to as pixels. Using MATLAB, this process becomes a 
bit easier as the principles involved in MATLAB for matrix 
manipulations can easily be used in image processing. One of 
the major applications of image processing is object or face 
recognition, in addition to which zooming in and zooming out 
an image is also a very important field of interest. In this 
paper, we are zooming in an image using cubic spline 
interpolation method. We programmed the method in 
MATLAB as m-file and tested it on 4 cameraman.tif . 

Cubic Spline Interpolation Method: This method relies on 
constructing a polynomial of low degree between each pair of 
known data points. If a first degree polynomial is used, it is 
called linear interpolation. For second and third degree 
polynomials, it is called quadratic and cubic splines 
respectively. The higher the degree of the spline, the smoother 
we get the curve. To obtain a smoother curve (or finer 
zooming), cubic splines are frequently recommended. The 
fundamental idea behind cubic spline interpolation is based 
drawing smoother curves through a number of points. This 
spline consists of weights attached to a flat surface at the 
points to be connected. A flexible strip is then bent across 
each of these weights, resulting in a pleasingly smooth curve. 
A cubic polynomial is the minimum order polynomial that 
guarantees the generation of C 0 , Cl, C2, curves. Therefore, it 
is the best known methods resort to piecewise cubic curve 
constructed by the individual third degree polynomials 
assigned to each subinterval. 

Mathematical Analysis of Cubic Spline Interpolation Process: 
The essential idea is to fit a piecewise function of the form 
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where Si(x) is a third degree polynomial defined by 
Si(x) = aix-Xjf + bi(x-Xif + Ci(x-Xi) + d t (l) 
x t are the data points, 
and i = 1,2,3, 

The fist and second derivatives of eq. (1) are given by 
Sj ’(x) = 3cii(x-Xi) 2 + 2b[(x-Xi) + Cj (2) 

Si”(x) = 6ai(x-Xi) + 2bj(3) 

where i being same as i = 1,2, n-1 

Our spline will need to conform to the following stipulations. 

1. The piecewise function Si(jc) will interpolate all data points. 

2. Si(v) will be continuous on the interval [xi,x n ] 

3. S[(x) will be continuous on the interval [x l9 x n ] 

4. S{\x) will be continuous on the interval [xl ,xn] 

Following the stipulations, we deduce the formulas for the 
coefficients as 
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■y ? i is the pixel color intensity at data point x- L 

After calculating the values of Mi using eq.(8), we calculate 
the values of a i? b i9 Ci and di for each pair of adjacent pixel and 
calculate the respective values for the zeroes filled in between 
them. The formula for filling the values in their respective 
positions was generated by hit and trial method. Now after 
contacting the image matrix and using the 4 imshow ’ 
command we get the zoomed image. 

For cubic spline, we calculate the values of Mi in terms of y 4 
using eq.(8) and after calculating the values M 2 , M 3 , . . . ,M n _i, 
we find the values of Ml and Mn using eq.(9) and eq. (10) 
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Results: Taking ‘cameraman.tif as sample grey scale image 
and the zooming intensity as 2, we can compare the results as 

the psnr value was found to be 148.48 and the execution 
time was 384.1070 seconds. 



Fig. 1 A part of cameraman.tif of size 50x50 



Fig. 2. 4x Zoomed in image of size 200x200 

So, we found that the zoomed image is much clear as 
compared to normally available zooming applications. 
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